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TECHNICAL MEMORANDUM 


LYMAN ALPHA SMM/UVSP ABSOLUTE CALIBRATION 
AND GEOCORONAL CORRECTION 

1. DEFINITIONS 

d. 

This report defines the instrument quantum efficiency K , where the index d 

s 

refers to detector number and s to slit number. This can be considered as the 
product of two functions 

K s =K A • 

where Q is a factor which is different for different slits. It has been assumed that 
s 

Q 2 h 1 ’ 


i.e., slit 2 as the reference slit. 

By definition, the function does not depend on slit but only on detector and 
time and will be expressed as (doy) where doy is a floating point number of days. 
Further K . and probably Q also depend on wavelength. 

Thus, K d can be considered as the ratio between detected photons, N, and the 
number of incident photons, Nj. The number of photons with frequencies between 

v - and v + , which enter the instrument during the gate time t g , is given by 


N. = t e a 

A V 


I dv 

_v 

hv 


> 


where h is the Planck constant, 1^ the radiation intensity, a the area of the observed 
solar surface, and w the solid angle which corresponds to the entrance of the instru- 
ment as seen from the Sun. 


The frequency can be taken out of the integral. The constant can be calculated 

2 

from the quotient of the telescope collector area (A = 66.4 cm [1]) and the square 
of the distance to the Sun (R 2 = 2.238E26 cm 2 ) to give 

o)= 2.967E-25 ster . 

The area of the observed solar surface can be expressed in square arc seconds 

s as 


/ 2 irR \ 2 _ 

a s ^ 360 x 3600 / 


s x 5.26E15 cm' 


so. 


_2 

oj a = s x 1.56E-9 cm ster 




and. 


aw 

he 


= 7.86E6 s . 


So for the number of incident photons, one writes (with X in cm) 


N i = t e s 7 * 86E6 x f ! v dv 
•'iv 


One may replace I 

•'A v 

equivalent. For Lya 


I dv with 

v 



1^ dX, provided the domains Av and AX are 


N. = t s 95.5 


I 


A v 


I dv 

v 


Since for slit 2 all lines (including Lya) are unresolved, i.e., much sharper 
than the slit, the maximum counts for that slit are proportional (linear response is 
assumed) to the integrated intensity in all cases. This is the main reason for using 
slit 2 as the reference slit. And since s = l”xl" , then 


2 


K 


d 



N 

max 


t e 7.86E6 X f I v dv 

Av 


y 


and for Lya 


v - ^max 

K d 

t 95.5 f I dv 

e J v 

A v 

When using other slits and unresolved lines, one can use similar formulae for 
but considering different s and X. 

When using other slits for resolved lines (not 1 or 2) , slits which have widths 
comparable to the line width, and for which the counts reflect the line profile (even 
when convolved), one can also calculate the integrated intensity by adding the counts 
at all scaning positions. But depending on the step size (n) and the effective slit 
width, one is counting many times the same photons (at the same frequencies) or 
missing some of them. In order to correct for this effect, one multiplies the counts 
by the multiplicity factor: 


„ _ n6 X 
m = 


AX 


eff 


where n6X is the wavelength displacement corresponding to each step. So one can 
write 


K = 


nS A 

t e S 7 • 86E6 X A 


En 


L 


I dv 

V 


For Lya one has 


R d _ n6 X . En 

s t e s 95.5 AX eff /* j dv 

J Av V 
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If one wants to obtain the absolute values of K;*, one can use the values for 

b 

the integrated intensities of the quiet Sun. See for example Vernazza and Reeves [2] 
which give fl^dv for many lines in the UV region of the spectrum. 


2. CALIBRATION 

Once the quantum efficiency has been determined, it can be used to calibrate, 
in absolute units, the profiles of resolved lines as well as the rasters. Assuming 
that 1^ is the intensity convolved with the slit response function ($ v ) normalized to 
1, one has: 


I c = 


L 


$ I dv or I , = 
V V A 


L 


'» : x dx 


with 


/ = / 

Av A A 


V X = 1 


Then, if one assumes that 1° does not vary over a n6 A interval (which is not 
always the case) , one finds for each step in wavelength 


N. = t s 7.86E6 A Av -- I 
l e ett v 


or 


N. = t s 7.86E6 A A A »» 1^ 
l e ett A 


where 


, v 4X eff e 4X eff 

4 v eff 1 X ' ~ • 

The coefficients C and C , are defined by the formula 

v A 


4 


where 


N = C I c = C, I® 

V V XX 


9 


C =K^ts 2.36E17 — ~ , 

v s e X 

and 

C .\ = K s t e s 7 - 86E6 x AA e ff ’ 

- 2 - 1 - 1 

(with X in cm and A X e ^ 1^ in erg cm s ster ). 

If using X in A and 1^ in erg cm 2 s ^ ster ^ A one has 

I = 3.33E-19 X 2 • I. 

v X 

For Lya 

I = 4.93E-13 I. 

v X 


3. SLIT PROFILE 


The only remaining problem is the definition of 4> and A X e ^ which can be 
obtained from the observed profile of a very thin line or a sharp edge, assuming one 
knows both integrated and peak intensity as well as the original line shape. The 
shape is not necessary if the original full width is much smaller than A X^. To 
obtain the function $ ^ , one can analyze the profile of a sharp unresolved line with 
n6 X << AX^. Under those conditions the observed profile is just the slit profile 

The slit effective width was taken from the table by Woodgate et al. [1]. The 
slit profiles seem to be well represented except for slits 1, 2, 18, 19, and 22 by the 
expression 


5 



$ 


X 


2 

irA X „„ 
eff 



A X \ 
AX eff / 


for 


0 


for 


AX 


< 


AX eff 

2 


AX 


AX 


eff 


For slits 1, 2, 18, 19, and 22, some theoretical diffraction curve may be used. 

The formula for the slit profiles can be used for deconvolution procedures, but 
one must remember that deconvolution is sometimes tricky and mathematically the 
solution of the problem is not unique. For this reason, it is more reliable to convolve 
the theoretical profile and then compare with the observed profile. 


4. THE GEOCORONAL CORRECTION FOR LYMAN ALPHA 

For the geocoronal correction of the Lyman alpha profiles, it is assumed that 
the intensity received by the instrument is related to the intensity , I® , from the Sun , 
by the formula 

I =1® exp ( - (J> t) 

v V * Y V 


where 


♦v =eXP (-(^) 2 ) ’ 

and <j> is the profile of the geocoronal absorption coefficient . t is the optical thick- 
ness of the geocorona at the center of the feature and A X Q is the Doppler width of 
the geocoronal absorption coefficient. For the conditions in the geocorona, the pro- 
file is well represented by a gaussian profile. 

This approach assumes that emissivity is negligible, which is true for solar 
observations, but is a poor assumption when observing very faint objects. Further, 
it assumes that the Doppler width is constant throughout the geocorona. 


6 


In the above formula, one has to consider also that both parameters t and A X^ 
can be functions of time, latitude, and longitude due to temporal and spatial changes 
in the geocorona. Also differential velocities can affect the present assumptions. 

The value of x to be used has to be considered a function of the zenith angle 
0 q, i.e., of the angle between the direction to the Sun and the Earth radius through 
the spacecraft position. 

The details of that function depend on the geometric model or the geocorona. 

It should be stressed that the dependence of x with Uq (v>q = cos6q) is not well 
represented by the formula 


T= T 0 /y 0 

because it can only be good for values of close to 1 and for a very thin geocorona, 
which is not the actual case. 

The model of the geocorona can be expressed as the function k (h) where k is 
the absorption coefficient at the center of the feature and h is the height (along the 
Earth radius). 

A model consisting of a constant value of k, which is zero except between 11 q 
and hp where hg is the spacecraft altitude, and h^ some other parametric value, 
was tried first. 

Since there are two experiments (V02345 and V02347) in which profiles have 
been obtained at various zenith angles, the authors have compared the variation of 
the feature depth with y Q and found that the model does not agree with the 
observations . 

Then a different model which seems to fit the observations was used. The 
model assumes an exponential variation of k, viz. 

/ r h - h o 

k = k 0 exp^-g- 

where Kq is the absorption coefficient at height 11 q above the Earth's surface and H 
is the scale height. This expression applies to the density in a hydrostatic, iso- 
thermal atmosphere. 
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A value of = 550 km was used in the calculations and it was assumed that 
SMM was in a circular orbit at that height. 

With dz being the element of distance from the SMM to the Sun, and y = cos 0 
the projection factor at height h for the ray with respect to the direction of the local 
radius, one has 


T 



oo 


z 0 


oo 



valid for 6 i 90 


Considering spherical geometry, 
siderations 


sin( ir- e Q ) 
h 


sin 6 



sin c 
z 


It follows that 


sin 



sin 0 q 


p is a function of h. From geometrical con- 
Sun > 



and that 



From the above expression, one sees that even when the formula of t seems to 
be the standard, y + jjq for all geometries except for the plane parallel case. 

From the above expressions one has 
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h - h o h 0 

or defining x = — g and 3 = -jj- 


u 

/ 

0 


t=k 0 HJ e' x (l + i) 


dx 


2 , x fo . x.-l 

“o + 6 (2 + e ) 


xTTT72 


If = 1, one finds 


K » H 

T o = T7 = K 0 H 


i.e., the optical thickness of the geocorona when observing along the Earth radius. 

No analytic expression for the integral has been found. It was calculated 

numerically and then tabulated as a function — (u n ,3). One interpolates in the 

T 0 U 

table for the different cases. 


The following values were obtained from the experimental data: 


A A D = 0.022 A , 
Tq = 0.63 , 


and 


3 = 2.33 


For the SMM orbit, some parameters of neutral hydrogen are calculated: 


a = -2L- -J— = 1.68E-18 
mc/iT D 


(“») 


= 9.28E-14 cm' 


T 0 - 2 
N„ = — = 6.78E12 cm 
ri a 
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H = 236 km 


T = 


_ n 2 

m H c 

2k 



1780 K° 


N H 

njj = -|j- = 2.87E5 cm 


m 


m 


kT 

m H gH 


6.36 


where a is the atomic cross section at line center, f the oscillator strength, N„ and 

II 

njj are the neutral hydrogen atoms column and particle densities, T the kinetic tem- 
perature (of neutral hydrogen atoms), the mass of hydrogen atoms, g the gravi- 
tational acceleration, c the velocity of light, k the Boltzmann constant, and the 
mean mass of particles with respect to hydrogen. 

From the above, one can correct the incident intensity for geocoronal absorption 
using the formula 

= I v exp(+<j>x) . 


The zenith angle can be calculated as suggested by Henze [ 3] , from the 
standard routines of Shine [4], knowing the time of the observation of the feature. 
Care has to be taken for long scans because the time delay during the experiment is 
sometimes important. 

The zenith angle cosine is given by the formula 
y q = cos 6 • cos a 

where 8 is the angle between the SMM orbital plane and the line of sight to the Sun 
and a is the angle around the orbit from the position of spacecraft orbital noon. 

a = 2ir(t - t n )/T 
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where T is the orbit period, t the time from orbital noon, and t the time of the 
observation. 

A detailed study of the time variations of the geocoronal parameters was not 
performed, but it was apparent for the 1980 data that these parameters did not change 
significantly. However, after the SMM repair and the observations began again in 
1984, the geocoronal parameters show larger variations, showing much deeper and 
perhaps even broader absorption and possibly a small turn to emission in the center 
of the Lyman alpha feature. 


5. CONVOLUTION OF GEOCORONAL FEATURE 

The above mentioned procedure for applying geocoronal correction is good only 
for slit 20 profiles, because in that case the slit width is smaller than the Doppler 
width of the feature. Since in the general case one has to correct some observed 
profile I®, one can write 

! v =/ = f <’ *0^ • 

A v A V 

If one applies the Fourier or Laplace transform to this expression, considering 
that the transform of the convolution of two functions is the product of their trans- 
form, one has 

Trans I® = Trans $ v [Trans exp(-<j> v x) + Trans I®] 


and 


Trans $ v • (Trans I®) = Trans I® - Trans [Trans exp ( - 1 }> ^ x) ] 

The correcting function is defined at a specific frequency as 

F = J $ v exp ( - <j> v t) dv , 

A v 
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CO 


and the corrected observed intensity I is 


I C0 = I c / F = 

V V 


'L 


$ I dv . 

V V 


The correcting function in turn can be obtained using the known slit profile 
and computing numerically the convolution. The complete expression is 


A A 


A + 


eff 


F = 


f 

2 2 r 

r/X - Y \1 

exp 

~ T 

. A X eff 

A 2 

A y COS 1 

" a * eff [ 

\ A X eff ). 

n - x o \ 2 

r p (^ D ) J 


dy 


where A Q is the wavelength at the center of the geocoronal feature and A is the 
wavelength at which one is evaluating the correction. 

From this formula, one sees that if AA 0ff is smaller than A A^ , the correction 
procedure from the previous section holds, i.e., F=exp(-<J> v x) . If this is not the 

A X eff 

case, one has to calculate the function F for all wavelengths between A Q ^ — and 

A A. 


A q + 


eff 


If the position of a ^ corresponds to one exact step , then one can write all the 
wavelengths in steps and get 


where 


g + 


‘'eff 


F = 


f 2 nn , 2 T 

w /H2_\ ‘ 

exp 

- T 

/ t ’‘eff [ 

« - ? 

( £ eff) 

1 y-so\ 2 

L exp U ) J 

A x e ff A X D 

A 

• g 0 

- X ° A 

o — 

' 6 A ’ 6 A ’ 

g 6 A 

- T x • 


dy 
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I 


Choosing g^ = 0 (i.e. , the zero wavelength in the geocoronal feature center) and 
defining 


p = .JL. = -Jr_ 

A eff A A eff 


and 


v = JL_= — A _ 
*eff A A eff 


one has 



v + 1/2 

2 2 

- cos [tt(v - p)] exp 

- 1/2 



dp 


The last integral is evaluated numerically for the integer values of g (steps in 
the wavelength drive) from - to + max(£ e ^/2, 2-e^). 


6. ABSOLUTE CALIBRATION OF DETECTOR 1 
AT LYMAN ALPHA 

The function has been evaluated by studying the rasters of the lambda max 
experiments with slit 2 and detector 1 at the Lyman alpha line. 

In those rasters, the experiments which show some quiet Sun close to the disk 
center have been selected. In the first observations it was usually impossible to find 
any image completely away from any active regions; in that case the images have been 
selected which show a portion of nonfacular areas and do not contain large 
filaments . 

Some tracings along diagonals and across the raster centers have been plotted 
for all the rasters from which the average counts have been taken for the whole 
raster and all the tracings. Comparing all average values one can easily recognize 
the quiet Sun areas which show constant average counts for all tracings, and the 
active region images which show higher values at some tracings. 

For the images containing active regions, average counts using selected portions 
of the rasters whose appearance was closer to the quiet Sun have been used. 
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The resulting values of quantum efficiency vary in steps with the time; that 
is, they are fairly constant for many days and then change suddenly in less than a 
few hours. All these steps (or degradation events), except the last case, seem to be 
correlated to the occurrence of large flares. It seems very likely that degradation 
could have occurred during the short time of a flare event. For large flares, the 
Lyman continuum intensity can rise by a factor of 5000 or more and the optical 
surfaces of the UVSP instrument may be adversely affected by this rapid and sig- 
nificant increase in UV radiation. 

The initial efficiency was (0) = 1.3E-3 which is close to the prelaunch 
estimate and is valid for the early experiments, up to doy = 61.11736. The next 
determination for doy = 67.98333 showed a new value (1) = 6.3E-4. The exact 
time of the change in sensitivity was traced using other slits and active region 
observations. The change occurs close to the flare event at doy = 66.18194, so 
t^ = 66.18194 is used. 

The next step occurred between doys 80.64651 and 85.64306 and can be traced 
to the compact, energetic flare observed at doy = 80.64651. The time is then t ^ = 
80.64651 and the efficiency after the event was (2) = 1.4E-4. 

Another step occurred on doy 89 and/or 90 in conjunction with the large flare 
observed in C IV line at 1548 A. The efficiency after this event was K^(3) = 5. IE- 5. 

The next change, and the largest one, took place between doys 97.4396 and 
97.5056, so t 4 = 97.50 was used. This event seems to be correlated with a large 
flare from Active Region 2373 at noon on April 6, 1980. The new efficiency was 
K x (4) = 1.7E-5. 

The last change in sensitivity occurred between the last data from 1980 and 
the first from 1984 and can be related to slow degradation. For the 1984 and following 
data, the efficiency K 1 (5) = 8.8E-6 is used. 

The authors looked at detector quantum efficiency for detectors 3 and 4 for the 
time interval 103.365 < t < 103.799 and obtained for this interval the following ratios: 

K 3 (t) / K x (4) = 0.79 ± 0.09 

K 4 (t) / K x (4) = 1.19 ± 0.06 . 
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Since there are insufficient observations, a detailed study of the time dependence of 

Kg and at the Lyman alpha wavelength was not performed. Data from other slits 

were used in order to obtain the coefficients Q . The best data were those from the 

s 

1984/1985 experiments, since they have quiet regions profiles taken with many slits. 

The values were also checked with earlier data to confirm the hypothesis of Q con- 

s 

stant with time. The values obtained for are: 

s 

Q 4 - 2.0 j Q-^2 = 1*0 > Q19 - 1.5 j Qgo - . 

For various reasons (including change in solar activity, location of slit relative to 
line center, and change in sensitivity), the uncertainty in the value of Q may be as 
large as a factor of 1.5. 

Another check on the consistency of the data comes from observations for 
235.4076 < t < 235.6430 using standard star Alpha Leonis on August 22, 1980 (Hayes 
and Shine [5]). For the different wavelength, the following was obtained: 

K^g(t) = 3.0E-4 at X = 1344 %. . 

7. FUTURE WORK 

A comprehensive collection of the calibration details of the UVSP instrument on 
the Solar Maximum Mission spacecraft is in progress. Additional slit, detector, and 
wavelength data will be tabulated and details concerning the polarization measurements 
and absolute calibration at other wavelengths will be represented in this follow-on 
work. In addition, the behavior of the Lyman alpha line in various solar features 
will also be described in a follow-on paper. The discussions contained in this report 
were preliminary to the papers which are to follow. 
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